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Abstract. We compute the Umit shapes of the Young diagrams of the minimal 
difference p partitions and provide a simple physical interpretation for the limit 
shapes. We also calculate the asymptotic distribution of the largest part of the 
Young diagram and show that the scaled distribution has a Gumbel form for all 
p. This Gumbel statistics for the largest part remains unchanged even for general 
partitions of the form E = nii^l" with v > Q where rii is the number of times 
the part i appears. 
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1. Introduction 

Exclusion statistics [U El O HI O [6l El [8] — a generalization of Bose and Fermi 
statistics — can be defined in the following thermodynamical sense. Let Z{f3, z) denote 
the grand partition function of a quantum gas of particles at inverse temperature (3 and 
fugacity z. Such a gas is said to obey exclusion statistics with parameter < p < 1, 
if ^(/3, z) can be expressed as an integral representation 



Jo 

where p(e) denotes a single particle density of states and the function yp{x), which 
encodes fractional statistics, is given by the solution of the equation 



In the cases p = and p = I, substituting yp{x) explicitly in ([T]) yield the standard 
grand partition functions of non-interacting bosons and fermions respectively. The 
fractional exclusion statistics with parameter < p < 1 (that corresponds to an 
interacting gas) smoothly interpolates between these two extreme cases. Two known 
microscopic quantum mechanical realizations of exclusion statistics are the Lowest 
Landau Level (LLL) anyon model [H [3] and the Calogero model [3 [7], with p(e) 
being, respectively, the LLL density of states and the free one-dimensional density of 
states. 

It is well known that a gas of non-interacting bosons {p — 0) or fermions 
{jp = 1) occupying a single particle equidistant spectrum both have a combinatorial 
interpretation in terms of the integer partition problem [S]. A partition of a positive 
integer E \s & decomposition of i? as a sum of a nonincreasing sequence of positive 
integers {/ij}, i.e., E = hj such that hj > /ij+i, for j = 1,2.... For example, 4 
can be partitioned in 5 ways: 4, 3 -I- 1, 2 -I- 2, 2 + 1 + 1, and 1 + 1-1-1 + 1. Partitions 
can be graphically represented by Young diagrams (also called Ferrers diagrams), 
where hj corresponds to the height of the j-th column (see figure [1]). In the Young 
diagram of a given partition of E, if ni denotes the number of columns having heights 
equal to i, then clearly E — UiSi — which can now be interpreted as the total 
energy of a non- interacting quantum gas of bosons where ei — i for i — 1, 2, . . . , oo 
represent equidistant single particle energy levels and Ui — 0, 1,2, . . . ,oo represents 
the occupation number of the i-th level (see figure [TJb)). On the other hand, if one 
expresses a positive integer E as a sum of strictly decreasing sequence of positive 
integers, i.e. E — J^j hj such that hj > hj+i (e.g. allowed partitions of 4 are: 4 and 
3 + 1), then the restricted partition problem corresponds to a non- interacting quantum 
gas of fermions, for which rii = 0, 1. In the partitioning problems if one restricts the 
number of summands to be iV, then clearly N = Ui represents the total number of 
particles. For example, if _E = 4 and N = 2, the allowed partitions are 3 + 1 and 2 + 2 
in the unrestricted problem, whereas the only allowed restricted partition is 3 + 1. The 
number p{E, N) of ways of partitioning E into N parts is simply the micro-canonical 
partition function of a gas of quantum particles with total energy E and total number 
of particles N: 



The grand partition functions, i.e., Z{P,z) = Sjv E_e -^^^ ^^p{E,N), for the 
unrestricted and restricted partitions are Z{P, z) = Y^iLii^ ~ ze~^'^)~^ and Z{(3, z) = 




(1) 



yp{x) -xyl 'P{x) = 1. 



(2) 




(3) 
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Figure 1. (a) All the Young diagrams for the partitions of 4. (b) The Young 
diagram of the partition 91 = 18 + 16 + 13 + 13 + 9 + 6 + 5 + 5 + 3 + 3, and 
the corresponding configuration of non-interacting bosons occupying energy levels 
Ei = i for i = 1,2, . . . , 18. 



ni^i(l + ^e-'^O and hence \nZ{(3, z) in the hmit /3 and p{e) — 1 reduce to ([T]) 
with p — Q and p — 1 respectively. 

Unhke Bose and Fermi statistics which describes non-interacting particles, for a 
quantum gas obeying exclusion statistics with parameter < p < 1, it is a priori 
not obvious how to provide a combinatorial description, since the underlying physical 
models with exclusion statistics describe interacting systems. However it has recently 
been shown [10] that a combinatorial description of exclusion statistics is possible in 
terms of a generalized partition problem known as the minimal difference p partition 
(MDP-p), which we will define in the next section. Even though the parameter p in 
MDP-p is an integer, in [TUj it has been shown that, when one analytically continues 
the results to non-integer values of p, for < p < 1, and in the limit f3 ^ 0, the 
MDP-p corresponds to a gas of quantum particles obeying exclusion statistics. This 
correspondence between exclusion statistics and MDP-p motivates us to investigate 
some other aspects of the MDP-p problem in this paper. 

2. Problems and outline 

In the MDP-p problem, a positive integer E is expressed as a sum of positive integers 
E = J2j^j such that hj — hj^i > p (see figure [2]). Therefore, p = corresponds 
to unrestricted partitions and p = 1 to restricted partitions into distinct parts. The 
shortest part in the MDP-p problem is usually taken to be > 1. However, for the 
calculation of certain specific quantities in this model, it is useful to consider a 
somewhat generalized version with the shortest part > s, where s is considered to 
be a variable. The grand partition function of this problem was obtained recently 
in [To], which is given by ([1]) with constant density of states p(e) = 1 and the lower 
limit of integration being s. 

One may also think of the MDP-p in terms of a quantum system consisting of 
equidistant energy levels = z for i = 1,2,...,cxd. Now a given height hj = i 
corresponds the energy level Ci — i and the number of columns with height i is the 
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Figure 2. A typical Young diagram for MDP-p problem. The thick solid border 
shows the height profile. Wh is the width of the Young diagram at a height h, 
i.e., Wfi is the number of columns whose heights > h. 



occupation number n^. Since the difference between two consecutive heights in the 
MDP-p must be at least p, the gap between two adjacent occupied energy levels must 
be at least p. Clearly for p = this gap is zero, and hence each level can be occupied 
by any number of particles (bosons). For p — I, each level can be occupied by at 
most one particle (fermions). Again for p > 1 a level can be occupied by at most 
one particle. However, in this case, when a energy level is occupied by a particle, the 
adjacent p ~ 1 levels must remain unoccupied. 

One major issue in the partition problem is to study the limit shape, i.e., the 
average height profile of an ensemble of Young diagrams with a fixed but large E. 
The shape (height profile) can be defined by the width Wh of the Young diagram at 
a height h (see figure [2]). In other words, Wh is the number of columns of the Young 
diagram whose height is greater than or equal to h. In this corresponding quantum 
system, Wh represents the total number of particles occupying energy levels above h. 

The height profile of the Young diagram of the unrestricted partition {p — 0) was 
first studied by Temperley, who was interested in determining the equilibrium profile 
of a simple cubic crystal grown from the corner of three walls at right angles. The 
two dimensional version of the problem — where walls (two) are along the horizontal 
and the vertical axes and E "bricks" (molecules) are packed into the first quadrant 
one by one such that each brick, when it is added, makes two contact along faces — 
corresponds to the p = partition problem. Temperley [TT] computed the equilibrium 
profile of this two dimensional crystal. More recently the investigation of the limit 
shape of random partitions has been developed extensively by Vershik [T^ |T31 |T3] 
and collaborators. The case of uniform random partitions was treated by Vershik 
who proved for the bosonic (p = 0) as well as the fermionic (p = 1) case that the 
rescaled hj^fE vs. Whl^fE curves converge to limiting curves when E — s- cx), and 
obtained these limit shapes explicitly. These results were extended by Romik [15) to 
the MDP-p for p — 1. In this paper we compute the following two quantities: 

(1) The limit shape of the Young diagrams of the MDP-p for any p, from which the 
previously obtained results for p = 0, 1, 2 follow as special cases. 
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(2) The distribution of the largest part of the Young diagrams of the MDP~p problem 
for all p, whereas the earlier result existed only for the p — case |16) . 

The average height profile (Wh) of the Young diagrams of the partitions of a given 
integer E is easier to compute in the grand canonical ensemble. Therefore one requires 
a restricted grand partition function Zh{l3, z) which counts the columns whose heights 
> h, and the full grand partition function Z{(3, z) which counts all the columns. From 
the restricted grand partition function one finds {Wh) — z-^ lnZh{(3, z)\z=i- For given 
large E, the parameter /? is fixed by the relation E — ~-§p In Z{(3, 1). 

On the other hand, to compute the number of partitions pp{E,l) of an integer 
E such that the largest part < Z, it is useful to consider the partition function 
= T.E'^'^'^PpiE'l) first. Formally Pp{E,l) can be obtained by inverting Zi{[3) 
with respect to (3, and for large E the asymptotic behavior of Pp{E, I) is obtained from 
the saddle point approximation, where the parameter P is fixed in terms of given E 
by the saddle point relation E — lnZ;(/3). 

Thus, it is useful to consider a more general restricted grand partition function 
Z(/3, z, Z, s) that counts the columns whose heights lie between s and /. All the other 
partition functions we need for our calculations can be obtained from Z{f3^ z, Z, s) by 
taking various limits on s and I. For example, by putting s — 1 and taking the limit 
I oo one obtains Z{j3,z). Similarly s — h and the limit / oo gives Zh{P,z) 
and putting s — 1 and z — 1 gives Zi{(3). As we will see later in ([T4| and ((23|) that 
/3 ~ E^^/^ for large E. Therefore, hereafter we will work in the limit /3 ^ 0. 

The rest of the paper is organized as follows. We first obtain the generalized 
grand partition function Z{(3,z,l,s) of the MDP-p problem in the next section. In 
sectionlHwe compute the limit shapes of the Young diagrams and also provide a simple 
physical interpretation of the result. In section [5] we calculate the distribution of the 
largest part of the MDP-p . Finally, we conclude with a summary and some remarks 
in section [6] 

3. Restricted grand partition function of MDP p problem 

Let Pp{E, N, I, s) be the number of ways of partitioning an integer E into N parts in 
the MDP-p problem such that the largest part is at most I and the smallest part is at 
least s, i.e., E = such that hi < I, /ij+i < hj — p for all j = 1, 2, . . . , — 1, 

and hpi > s. Then clearly, [pp{E,N,l,s) — pp{E,N,l — l,s)] gives the number of 
MDP-p of E, such that the largest part is exactly equal to I, and smallest part is at 
least s. Now, by eliminating the first part hi = I from the partition one immediately 
realizes that the above number is precisely Pp{E — I, N — 1,1 — p, s), i.e., the number 
of MDP-p oi E — I into A^ — 1 parts such that the largest part is at most / — p and 
the smallest part is at least s. Therefore, one has the recursion relation 

Pp{E, N, I, s) = Pp{E, N,l-l,s)+Pp{E~l,N^l,l-p, s). (4) 

Following similar reasoning one can also derive another recursion relation in terms of 
the smallest part s, 

Pp{E, N, I, s) = pp{E, N,l,s + 1) + pp{E ~s,N~l,l,s+ p). (5) 

It follows from ^ and ([5]) that the grand partition function Z{f3,z,l,s) ~ 
Sat J2e Ppi^j Af, Z, s) satisfies the recursion relations: 

Z(/3, z, I, s) = Z(/3, z,l-l, s) + ze->^^Z{(3, z,l~p, s), (6) 

Z{P, z, I, s) = Z{P, z,l,s+l) + ze-^'Z{P, z,l,s+p). (7) 
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From these equations, it is evident that in the scahng hmit /3 — s- 0, and both s and / 
large, the correct scaling variables are ps and (31, so that /3s and (31 remain finite. One 
knows from the statistical mechanics that the free energy Z{(3, z,l, s) becomes a 
function of the only the scaling variables in the limit P 0. Therefore in this limit it 
is natural to expect 

Z(/3,z,/,s) « exp ( i$(/3;,^s,z) ) . (8) 



Now to determine the scaling function (f>{(3l, Ps, z), we substitute the ansatz ^ in 
([U and (O, and then expand ^{Pl — P,Ps,z) and ^{Pl — Pp,Ps,z) about pi, and 
^{Pl, Ps + P, z) and <I>(/3/, Ps + Pp, z) about /3s, respectively in Taylor series up to first 
order, which yields the equations: 

d 

exp(-$/3z) + ze"^'exp(-p$0,) = 1> where ^f}i = —<^{u,(3s,z) (9) 



u=f3l 



d 

exp($/3s) +ze-'^"exp(p$^,) =1, where = — d, z) . (10) 

av v=f3s 

It is evident from Q and pO]) . that <i>^/ and $/3s are function of the arguments ze~'^' 
and ze"''" respectively, and the solutions are 

$/3; In Vp (ze"''') and $;3s = - hi yp {ze''^') (11) 

where yp(x) satisfies the equation yp{x) — xyp~P{x) = 1, which is the same equation 
([2]) one encounters in exclusion statistics. Equation (fTT|) implies, $(m, u, z) = 
J" Inyp {ze~'^) de. Therefore, H]) yields 

1 r/^' 

\nZ{P,z,l,s)^- \nyp{ze-') de, (12) 

i.e. ((!]) with constant density of states p{e) — 1, and the lower and upper limits 
of integration being s and I respectively. This is the key equation, using which we 
compute the limit shapes and the largest parts of the Young diagrams in section 2] 
and section [5] respectively. The limit (31 ^ oo also provides a simpler derivation of 
an earlier result [lOj . which showed a link between the exclusion statistics and the 
MDP-p problem. 



4. Limit shapes of Young diagrams 

Let us consider all the MDP-p of an integer E with uniform measure. Then the 
number of columns having height between s and I, averaged over all the Young 
diagrams of the MDP-p of E, is obtained from as 

(iVi(z)) =z^lnZ(/3,z,Z,s) = ^ [iny^ze-'^^) - Iny^ (zg-^')] • (13) 

Now to obtain the parameter (3 in terms of the given large integer E one again uses 
(|12p with the limits pi — > oo, /3s 0, and z ~ 1, i.e., 

ii; = -^lnZ(/3, 1,00,0) = ^, where b^p) = J^°" Iny^ {e-') de (14) 

is a constant which depends on the parameter p. 
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Figure 3. Limit shapes for the minimal difference p partitions with p = 0, 1,2, 
and 3, where 6(0) = tt/VG, b{l) = tt/Vu, 6(2) = n/VTE, and 6(3) = 0.752617 . . .. 



P{Wh) =lnyp{e-P'') , where p - 



(15) 



The average shape or the height profile of the Young diagrams (Wh) is simply 
given by (|13p with s = h, I ^ oo and z = 1, i.e., 

For instance for p = 0, 1 and 2, solving ([2]) yields ynix) — 1/(1 — a;), yi{x) — (1+x), and 
2/2(3;) = [1 + V 1 + 4xJ /2 respectively. From which using (fT4)) one finds 6(0) — tt/\/6, 
6(1) = tt/^/12 and 6(2) = tt/vTs in agreement with the earlier known results I15j . 
The fiuctuation about the average shape can be computed from (fT2t using 



{W,,f ^z — z—lnZ{(3,z,^,h) 



dz dz 



which gives 



-l3h 



(e-'5^) 



(16) 



(17) 



where y'p{x) denotes the derivative of yp{x) with respect to its argument. This 
formula shows that the random variable fiWh is strongly peaked around its mean 
value. Therefore, the curve Wh/^/E as a function of h/\/E converges to a limit curve 
when E 00 (strictly speaking, to prove the existence of a limit curve, one needs to 
show that all the moments around the mean vanish when E —^ 00, which Vershik |13j 
showed for p = and p = I). Therefore hereafter we may replace {f3Wh) by f3Wh- 
Using (21) and (fT5)) . one can express h in terms of Wh as, 

1- - 'PW,\_pWf^_ (18) 



■In 1 



Introducing the scaling variables x — Wh/VE and y — h/^/E, using (fT4)) and taking 
E' — !■ 00, yields the equation of the limit shape 



y ■ 



-xb{p) 



) -px. 



(19) 



Figure [3] shows the limit shapes for the MDP-p with p = 0, 1, 2, and 3. 
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(a) 



(b) 



(c) 



(d) 



Figure 4. (a) Transposed Young diagram for tlie unrestricted partition 91 = 
18 + 16 + 13 + 13 + 9 + 6 + 5 + 5 + 3 + 3. (b) Represents (a) in terms of 
non-interacting bosons (represented by O ) occupying energy levels ti = i for 
i = 1,2,..., 10. (c) The configuration obtained from the bosonic configuration 
(b) by transferring particles from the higher levels to the lower levels such that 
in the final configuration all the levels below the highest occupied level £5 = 5 
receive 2 new particles (represented by • ) each, where O represents the particles 
originally present in the initial bosonic configuration, (d) The Young diagram 
corresponding to the configuration (c). This is the transposed Young diagram 
of the partition 49 = 18 + 14 + 9 + 7 + 1, in the minimal difference 2 partition 
problem. 



Equation (jlSp has a simple physical interpretation which we explain below. For 
p = 0, any transposed Young diagram (see figure [3]) provides a valid unrestricted 
partition. Therefore the transposed diagram also corresponds to a non interacting 
system of bosons occupying single particle equidistant energy levels. However this is 
no more true when p > 0. In this case the transposed Young diagram(see figure |4]) 
corresponds to a quantum system where there is a certain energy level (which differs 
from one realization to another) which is occupied by at least one particle, and above 
which all the levels are empty, and below which each of the levels must be occupied 
by at least p particles. Therefore, h in the limit shape expression p8)) represents the 
number of particles above the energy level Wh- For bosons with total energy this 
number is precisely given by (jlSp with p = and (3 has to be determined in terms 
of E. Now, a configuration for p > 0, can be obtained from a bosonic configuration 
by transferring particles from the higher energy levels to the lower ones such that, 
in the final configuration, levels below the highest occupied level (which has at least 
one particle) receive exactly p new particles each. Clearly, in the final configuration 
obtained by this procedure, each of the levels below the highest occupied level has at 
least p particles. However, since transferring a particle from a higher energy level to 
a lower one decreases energy of the system, to obtain a configuration for p > with 
energy E requires the initial bosonic configuration to be at a higher energy (i.e., lower 
inverse temperature (3) than E. Now, while going from a initial bosonic configuration 
to a configuration for p > 0, one transfers total oipWh particles from the levels above 
Wh to below (i.e., p particles to each level), the average number number of particles 
above level Wh decreases from the corresponding bosonic system (p — 0) precisely 
by pWh, which is exactly the content of (HH]) . In fact, (3 in (|18p can directly be 
determined by using condition h > and the normalization J^'^ h{Wh) dWh = E, 
where Wj^ is the solution of the equation h{Wj^) — 0. Writing exp{l3W^) — y* , it 
satisfies y* — y*^~^ = 1, and in terms of y* one finds 

P=M with b'ip) = ^-U,il/yn-^{lnyr, (20) 
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where Li2(z) = X^fc^i ^'^^ ^ ^he dilogarithm function. The expression for b{p) in 
(|20p also follows directly from its integral representation given in (fT5|) . 

5. Largest part of Young diagrams 

Equation also allows one to compute the distribution of the largest part (i.e., 
the largest height in the Young diagram) in the MDP-p problem. Let Pp{E, I) be the 
number of partitions of the integer E in MDP-p problem, such that the largest part is 
at most I. Clearly, Pp{E) = Pp{E, I — > oo) gives the total number of partitions of E and 
since the partitions are distributed with a uniform measure Cp{l\E) = Pp{E, I)/ pp{E) 
gives the cumulative distribution of the largest height /. Note that the partition 
function Z/(/3) = e~^^pp{E, I) in the limit /3 ^ is obtained from by simply 
taking the limit f3s —^ and z = 1. Therefore, formally inverting the Laplace transform 
(in the limit /3 — )■ 0, the sum over E in the partition function of Zi{j3) can be replaced 
by an integral), one can write 

Pp(E,l) = — exp[5£;,(/9)] d/3, (21) 

'^'^1 J~j—ioa 

where 7 is a real constant chosen such that all singularities of integrand are to the left 
of the vertical contour in the complex plane, and the action 

1 

SeAP) =PE+- \nyp (e"^) de. (22) 
P Jo 

For large E, the leading asymptotic behavior of pp{E,l) can be obtained from the 
saddle point approximation. Maximizing the action with respect to f3, i.e., setting 
dS/d(3 = gives the saddle point equation 

(3^E^ \nyp (e"^) de - pi In yp {e-^') . (23) 
Jo 

For large E, the saddle point (3* is obtained implicitly solving the above equation and 
by substituting it back in the action Se,i{P*)- Thus, to leading order, 

Pp{E,l)^e^p[SE,i{P*)], (24) 
where Se,i{P*) can be written as 

lnj/p(e-^) de-/3*nnyp(e-^ . (25) 



It is evident from the above equations that, in terms of I and E, one has the scaling 
form SeAP*) = ^J~E gp{l I \fE\ where the scaling function gp(x) can be determined as 
follows. We set Ij^fE = x and fi*l = Hp{x). In terms of these scaling variables, from 
the saddle point solution of ([23ll and the entropy ([25| one has 

H'^(x) fHpix:) 

-j^ = Int/p (e-^) de - Hp{x) \nyp [e-""^^^) , (26) 

and gp{x) = 2^^^^ + xlny^ (e-^-^^)) , (27) 

respectively. Thus, given x, one has to find Hp{x) by implicitly solving p6p . then 
substitute it back in (|27p to get gp{x), and finally 



PpiE, I) w exp 



(28) 
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For large x, using 

gp{x) sa 2b{p) - 



and (P7)l . it can be shown that 
1 



Hp) 



exp [—b{p)x] 



as X 



oo, 



(29) 



where b{p) is given in HH) and (gOl)- Thus, from ^E^, Pp{E) = Pp{E,l ^ oo) ^ 
eyip[2b{p)\/E] to leading order for large E, which is the generalization of the Hardy- 
Ramanujan formula [17] for po{E), provided by Meinardus [18]. The normalized 
cumulative distribution of I, i.e., Cp{l\E) = pp{E,l)/pp(E), for large E and I ^ VE, 
is therefore 



Cp{l\E) « exp 



E 



exp 



Hp) 



= F 



E 



l-r{E) , (30) 



where the characteristic value of I is r{E) — [\/E/b{p)\ \n{\/ E / b{p)) , and the scaling 
function has the Gumbel form, F{z) = exp[— exp[— z]]. The result for the p = case, 
i.e., for Cq{1\E), was first derived Erdos and Lehner [IB]. Equation provides a 
generalization of their result, which is valid for all p. The probability distribution 
Pp{l\E) = Cp{l\E) - Cp{l - l\E) « dCp{l\E)/dl, obtained from (fSO]), 



Pp{l\E) 



Hp) 



E 



F' 



fe(p) 
E 



I — 1*{E) ) , where F'{z) — exp[— z — cxp[— 



(31) 



is highly asymmetric around the peak at I = 1*{E). This limiting distribution describes 
the probability of typical fluctuations of 0{-\/E) of the random variable / around the 
peak 1*{E). 



6. Summary and remarks 



In summary, we have obtained a generalized grand partition function for the minimal 
difference p partition (MDP~p) of a positive integer E, where smallest part is at least 
s and largest part is at most I, in the scaling limit /3 (x E~^^^ 0, in terms of the 
scaling variables f3l and /3s. The limit f3l —f oo also provides a simpler derivation of 
an earlier result [10], which showed a link between the exclusion statistics and the 
MDP-p problem, by showing that both problems are described by the same grand 
partition function in the limit f3 0. Using the grand partition function we have 
computed the limiting shape of the Young diagram of the MDP-p problem for all p, 
and also provided a simple physical interpretation of the result. Although the Young 
diagram is defined only for integer values of p, one can analytically continue the 
expression ([TS]) for the width Wh of the Young diagram to noninteger values of p. For 
< p < 1, Wh corresponds to the number of particles each of which has energy at least 
h, in a system where the particles obeys exclusion statistics. We have also obtained 
the asymptotic distribution of the largest part of the Young diagram and showed that 
the scaled distribution has a Gumbel form for all p. When one analytically continues, 
for < p < 1, the largest part corresponds to the highest occupied energy level in 
exclusion statistics. 

Note that for p = 0, the transposed Young diagram of a given partition gives 
another valid p = partition. This symmetry implies that the statistics of the 
largest part is the same as the statistics of the number of parts in the p = partition 
problem. The distribution of the number of parts for p = was computed by Erdos 
and Lehner [16] and in the appropriate scaling limit it has a Gumbel form. However, 
the symmetry between the number of parts and the largest part no longer holds when 
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p > 0, where the distribution of the number of parts become Gaussian (see [lOj and 
references therein). 

Recently, the statistics of the number of parts for a general partitions of the form 
E = Uii^/'^ that corresponds to having a power-law density of states, p{t) ~ e'^^^, 
has been studied [19] in the bosonic sector (p = 0). Clearly, v = 1 corresponds to 
the usual unrestricted partition problem, where the number of parts obey Gumbel 
statistics. Interestingly, for 7^ 1, the authors in [19] also obtained the other two 
universal distribution laws of extreme value statistics, namely the Frechet and WeibuU 
distributions for < < 1 and v > 1 respectively. 

Therefore, the general partition problem can be defined in the parameter space 
of (i', p) with ly > and p > 0. In this parameter space the point {i^ = l,p = 0) is 
a very special one at which both the number of parts and the largest part obey the 
same statistics given by the Gumbel distribution. Along the line v = 1, the limiting 
distribution of the number of parts becomes Gaussian as soon as p > 0, whereas the 
limiting distribution of the largest part remains Gumbel for all p, as we have shown 
in this paper. On the other hand, along the p = line, for the number of parts one 
finds [H] all the three universal laws of the extreme value statistics, for the parameter 
< V < I, v = 1, and v > 1 respectively. Therefore, it is interesting to ask whether 
there is any region in the (z^, p) parameter space, where the largest part obeys another 
statistics than the Gumbel one. The answer is negative. For a general density of 
states, (|22p includes a factor of /5(e//3) in the integrand. Following the similar steps 
provided afterwards, it can be shown that even for the power-law density of states 
p(e) ^ e'^^^, the scaled distribution of the largest part remains Gumbel in the whole 
{v,p) plane. Thus, the largest part obeys a more robust law, in contrast to the number 
of parts. 



Note added in proof 



We thank K. Hikami for pointing out [20] in which the author obtained the solution 
of a recursion relation similar to ([6]) with s = for arbitrary /3. However, for the 
purpose of this paper we require the solution only in the limit /3 — > 0. In this limit it 
is simpler to obtain it using the method presented in this paper rather than obtaining 
by taking the limit /3 ^ in the solution of [5D]. The average occupation number at 
a level i for the exclusion statistics has been studied in [51|31 IH [HI HI] , which also can 
be obtained from simply through 



which via elementary algebra yields 

1 



d_ 

di 



{Ni{z)) 



l=i 



(32) 



(33) 
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